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We present a new covariant and gauge-invariant perturbation formalism for dealing with space- 
times having spherical symmetry (or some preferred spatial direction) in the background, and apply 
it to the case of gravitational wave propagation in a Schwarzschild black hole spacetime. The 1+3 
covariant approach is extended to a '1+1+2 covariant sheet' formalism by introducing a radial unit 
vector in addition to the timelike congruence, and decomposing all covariant quantities with respect 
to this. The background Schwarzschild solution is discussed and a covariant characterisation is 
given. We give the full first-order system of linearised 1+1+2 covariant equations, and we show 
how, by introducing (time and spherical) harmonic functions, these may be reduced to a system of 
first-order ordinary differential equations and algebraic constraints for the 1+1+2 variables which 
may be solved straightforwardly. We show how both the odd and even parity perturbations may 
be unified by the discovery of a covariant, frame- and gauge-invariant, transverse-traceless tensor 
describing gravitational waves, which satisfies a covariant wave equation equivalent to the Regge- 
Wheeler equation for both even and odd parity perturbations. We show how the Zerilli equation 
may be derived from this tensor, and derive a similar transverse traceless tensor equivalent to this 
equation. The so-called 'special' quasinormal modes with purely imaginary frequency emerge natu- 
rally. The significance of the degrees of freedom in the choice of the two frame vectors is discussed, 
and we demonstrate that, for a certain frame choice, the underlying dynamics is governed purely 
by the Regge- Wheeler tensor. The two transverse-traceless Weyl tensors which carry the curvature 
of gravitational waves are discussed, and we give the closed system of four first-order ordinary dif- 
ferential equations describing their propagation. Finally, we consider the extension of this work to 
the study of gravitational waves in other astrophysical situations. 



I. INTRODUCTION 



The 1+3 covariant approach has proven to be an extremely useful technique for developing a detailed understanding 
of many aspects of relativistic cosmology, both in terms of fully nonlinear GR effects and through the application 
of the gauge- invariant, covariant perturbation formalism (see for example) to the formation and evolution of 
density perturbations Q in the universe and to the physics of the cosmic microwave background 0, Q , amongst other 
things (see pj for a comprehensive review). Rs strength in cosmological applications lies in the fact that it is well 
adapted to the system it is describing: all essential information can be captured in a set of (1+3) covariant variables 
(defined with respect to a preferred timelike observer congruence) , that have an immediate physical and geometrical 
significance. They satisfy a set of evolution and constraint equations, derived from Einstein's field equations, and 
the Bianchi and Ricci identities, which form a closed system of equations when an equation of state for the matter 
is chosen. The covariant and gauge- invariant linearisation procedure is easy and transparent: it consists of deciding 
which variables are 'first order' (or 'of order e') and those which are 'zeroth order' - i.e., those which do not vanish 
in the background, which is usually a Friedmann-Lemaitre-Robertson- Walker (FLRW) model. Products of first-order 
quantities can then be ignored in the equations. Whenever the background is a homogeneous and isotropic FLRW 
model, all projected vectors and tensors are first-order, so there is no vector-tensor and tensor-tensor coupling in 
the equations. Harmonic functions can then be introduced which re-write the equations in scalar form; the resulting 
system is then in the form of algebraic constraints and some first-order ordinary differential equations; the solution 
is then straightforward. The key point of the approach is that it deals with physically or geometrically relevant 
quantities, such as the fractional density gradient, V a , or the electric and magnetic parts of the Weyl tensor, E a j, and 
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H a b, respectively, which represent the non-local parts of the gravitational field, and which describe, amongst other 
things, the propagation of gravitational waves (GW). The variables are also (coordinate) gauge invariant (although 
there is a frame-gauge freedom in the choice of u a - see below). 

The aim of this paper is to extend the gauge-invariant, covariant perturbation formalism to an astrophysical setting. 
The 1+3 approach is not appropriate for many situations where such techniques would seem highly desirable: when 
the spacetime in inhomogeneous, for example, the 1+3 equations usually become intractable. However, by introducing 
an additional frame vector, assuming that the background spacetime has some preferred spatial direction (such as in 
spherically symmetric, or more general locally rotationally symmetric spacetimes, or G2 spacetimes) we can in many 
cases recover all of the advantages of the 1+3 equations, but in a 1+1+2 covariant framework 5j. In this paper we 
introduce the 1+1+2 formalism, and apply it to linear perturbations of a Schwarzschild spacetime. Not only is this 
a first step in applying the procedure to more general astrophysical situations (such as perturbations of the interior 
of compact objects, collapsing and exploding stars, etc.), it also represents an important field of study in itself: with 
the development of large gravitational wave detectors (e.g., 0) an improved understanding of the problem of GW 
propagation around compact objects is certainly timely. The power of the 1+1+2 technique is clearly shown by the 
significant results we are able to obtain, relatively simply. For example, we show here that the full description of 
gravitational waves around a Schwarzschild black hole is governed by closed covariant wave equation, unifying both 
parities in a single covariantly defined gauge- and frame-invariant transverse-traceless 2-tensor, W a b- 

Linear perturbations of Schwarzschild black holes have been conventionally studied through perturbations of the 
metric tensor (or via the Newman-Penrose formalism 7]). In the metric approach, fluctuations of the spacetime 
geometry are characterised by perturbations in the metric tensor; these fluctuations are determined by closed wave 
equations - the Regge- Wheeler equation for odd parity perturbations Q and the Zerilli equation in the even parity 
regime t 9j. These wave equations act on linear combinations of the functions (and their derivatives) appearing in 
the perturbed metric, but these functions do not determine directly the gravitational waves which they represent; 
a general coordinate transformation would preserve neither the perturbation functions themselves, nor the wave 
equations which they satisfy. The approach we develop here is completely covariant, so such issues do not arise. 
Instead, corresponding wave equations we derive here are formed from covariant and gauge-invariant variables which 
have a physical significance; furthermore, they do not require a harmonic splitting for their derivation. 

The formalism we develop here relies on a further splitting of the spacetime using a radial vector n a , in addition to 
the usual splitting with the timelike vector u a used in the 1+3 approach. We split the Ricci and Bianchi identities 
using u a and n a into a coupled set of first order differential equations, plus some constraints. The differential operators 
we use are along the two vector fields which give us a set of evolution and propagation (along the 'radial' direction) 
equations, while a derivative formed from a projection orthogonal to u a and n a gives a small number of constraints. 
The differential equations involve the covariant variables derived from splitting the Weyl tensor (and more generally 
the Ricci tensor, but this is zero here as we only consider vacuum perturbations), and the kinematics of u a and n a . As 
our background is static and spherically symmetric, we may use harmonic functions for our evolution and projected 
derivatives, putting our equations into the form of a first order system of ordinary DE's and constraints, which can 
then be tackled relatively easily. 

Previously, 0, 0, Il2t fli^ (and references therein) developed approaches to stellar and black hole perturbations 
similar to the method presented here in the sense that they use two orthogonal vectors to form their time and space 
derivatives. These approaches are fundamentally different from our approach, however, in that they formulate their 
differential equations as second order PDE's derived from Einstein's field equations (EFE), the solutions of which give 
the metric functions (as in all metric perturbation approaches). On the other hand, our system of DE's is manifestly 
first order, as it is derived from the first order Ricci and Bianchi identities 1 , as it involves physical or geometric 
quantities, and not the metric functions. Second order wave equations may be derived if desired. This is one of the 
key properties of covariant or tetrad methods. This change in derivative level is conceptually analogous to the change 
in going from the Lagrangian to Hamiltonian formulations of classical mechanics, from configuration space to phase 
space. 

The layout of the paper is as follows: in the following section we discuss the merits of a 1+1+2 decomposition 
of the field equations and set out the 1+1+2 covariant formalism. Then, in section ITTT1 we present the full set of 
1+1+2 covariant, gauge- invariant, first-order equations, linearised about a Schwarzschild background and introduce 
the (spherical and temporal) harmonics on the 'sheet', which enable the equations to be reduced to a set of coupled 
ODEs for the 1+1+2 covariant variables. In section llVl we prove the existence of a transverse-traceless (TT) tensor 
that satisfies a closed wave equation equivalent to the Regge- Wheeler equation, valid for harmonics of either parity. 



We use the Ricci identities applied to our vectors u a and n a , so they are second order DE's of u a and n a ; however, our covariant objects 
are different projections of V a Uj, and Va^ij, so that the Ricci identities are first order DE's in these covariant objects. 
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Following this we discuss the even parity variable which satisfies a wave equation equivalent to the Zerilli equation; 
we demonstrate here the existence of an odd parity counterpart, but defer the derivation of the 'Zerilli tensor' until 
later, Sec. IV Cl Then, in section we describe in detail the (matrix-based) method of solution of the linear, first- 
order system of ODEs for the harmonic components, emphasising the significance of the freedom to choose the frame 
vectors and showing that with an appropriate choice we can reduce the whole solution for both parities to a single 
2-dimensional ODE. The wave equations for the TT electric and magnetic Weyl tensors are also presented, and the 
closed four-dimensional ODE which they also satisfy is discussed. Finally, we discuss the results we have obtained. 
We follow the notation and conventions of fil. 



II. THE 1+1+2 COVARIANT SHEET APPROACH 

Before setting out the principles and equations of the 1+1+2 covariant formalism it will be illuminating to examine 
the 1+3 approach to see where its strengths and weaknesses lie, and why it is so successful in a cosmological setting, 
but is less useful in other situations, such as in the study of gravitational radiation in a Schwarzschild spacetime which 
we consider here. In the process we hope to indicate that the 1+1+2 formalism neatly fills a gap between the 1+3 
and tetrad approaches. 

A. 1+3 covariant perturbation theory: why it works in cosmology, but not for black holes 

In a nutshell, the 1+3 formalism is successful in cosmological applications because the assumed spatial homogeneity 
and isotropy of the background spacetime means that the only essential coordinate is time: the introduction of an 
appropriate timclikc observer congruence u a allows the full structure of the spacetime to be described solely in terms 
of ordinary differential equations involving (1+3) scalar quantities because all spatial derivatives, spatial projections 
of vectors, and projected, symmetric, trace-free parts of tensors must vanish by symmetry. 

On the other hand, when the spacetime is not homogeneous and isotropic the resulting equations are not simple 
ODEs. Consider, for example, a family of static observers around a Schwarzschild black hole, that is, observers on the 
congruence u a parallel to the timelike conformal Killing vector. Then u a has zero rotation, shear and expansion, but 
has non-zero acceleration ii a (reflecting the fact that a force must be applied to prevent infall) , and the electric Weyl 
curvature, E ab , is non-zero, while all other covariant quantities are zero. Relative to these observers, E ab measures 
non-local gravitational effects: in this case the (time-independent) radial tidal forces only, which can be described by 
a single function of distance from the hole - thanks to the spherical symmetry we can think of it is a tensor describing 
an essentially scalar phenomenon. The 1+3 covariant equations describing the spacetime are 2 



E {ab) = 0, (1) 

div?i + ?i 2 = 0, (2) 

curl w a = 0, (3) 

Eab - D( Q U 6 ) - U( a U b ) = 0, (4) 

cur\E ab + 2e cd{a u c E b) d = 0, (5) 

div£ a = 0. (6) 



This is already a formidable set of tensor equations; indeed, their solution is basically an intractable problem unless 
we introduce a full tetrad or revert to a metric based approach. These problems become even more severe when we 
consider the perturbed spacetime. 

In the 1+3 approach this is achieved by assuming that in general all of the 1+3 covariant quantities are non-zero, 
but that any that vanish in the background are small - small enough that we can neglect products of such terms. To 
get an idea of the horrendous nature of the equations, consider the (gauge-dependent) wave equation for E ab , which 
gives information about gravitational waves (or linear dynamical tidal forces) : 

E (ab) - D 2 E ab = -D c D (a E b)c + 5u c D c E ab - 2u c D {a E b)c - 3E c(a D c u b) - E c(a D b) u c 



2 We use the standard notation whereby a dot represents differentiation along the observers' four-velocity u a , ip a --b = " c V c i/>a---bi and 
D a is a derivative in the rest space of the observers; D c i/> a -.f> = h^h^ ■ ■ ■ (i/Vrf,-/) where = g a ^ + u a ui, is the usual projection 
tensor orthogonal to u a . We use angled brackets on indices to donate the projected, symmetric and trace-free part of a tensor. 
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+ \e cd{a cur\ ii c E b) d + 2{d\v ii + u 2 ) E ab + 3E c {a E b)c - 6u c u (a E b)c . (7) 

This also contains information about non-linear tidal forces, which don't propagate at the speed of light, by virtue 
of the presence of the parts of E ab that do not vanish in the background. But there is no way to separate the two 
physical effects. In addition, how could we solve this wave equation? In contrast to FLRW models, which have only 
scalars describing the model after a 1+3 decomposition, it is the presence of non-zero vectors and tensors in the 
background spacetime which makes a black hole impossible to deal with in the 1+3 approach: all the equations have 
vector-tensor and tensor-tensor coupling in them, rendering them intractable. 

The key, then, for the covariant perturbation approach lies in the fact that, in the background, the congruence u a is 
orthogonal to 3-surfaces of homogeneity and isotropy. In the case of inhomogeneous, spherically symmetric systems, 
projection via u a is simply not enough. Another vector field is required that is orthogonal to homogeneous and 
isotropic surfaces. Clearly, after an appropriate projection with u a , such surfaces are provided by spheres surrounding 
the centre of symmetry, and the vector orthogonal to these is a radial vector. We turn now to developing such a 
formalism. 



B. The 1+1+2 formalism 

In the 1+3 approach, a timelike threading vector field u a (u a u a = —1) is introduced, representing the observers' 
congruence. Given this vector field, the projection tensor h b = g b + u a u b is introduced, which projects all vectors 
and tensors orthogonal to u a . Using h ab , any 4-vector may be split into a (1+3 scalar) part parallel to u a and a 
(3- vector) part orthogonal to u a . Any second rank tensor may be covariantly and irreducibly split into scalar, vector 
and projected, symmetric, trace-free (PSTF) 3-tensor parts, which requires the alternating tensor e abc = ^rjdq^jM- 
Tensors of higher rank may be similarly split, but are rarely used (an important exception being CMB physics |3j,|j|). 
These are the fundamental quantities describing the spacetime, after the introduction of u a . 

We now introduce another vector field and perform another split, but this time of the 1+3 equations. The '1+1+2' 
decomposition we develop here has been partially studied before, mostly in the context of symmetries of solutions 
of the EFE dEHTil. It was introduced by flfij ] and further developed in [HI Il7j|. However, there are importances 
differences with the work presented here. In the following we assume the 1+3 covariant split of the equations (as 
given in 0, for example), with all tensors split into scalars, vectors and PSTF tensors with respect to u a . 

Take a unit vector n a orthogonal to u a : n a n a — 1, u a n a = 0, and define the projection tensor 

N b EE h a b - n a n b = g a b + u a u b - n a n b , (8) 

which projects vectors orthogonal to n a (and u a ): n a N ab = = u a N ab , onto 2-surfaces (N a a = 2) which we refer to 
as the 'sheet' (to carry the sewing analogy of the threading approach into the realm of the ridiculous). 

Any 3-vector ip a can now be irreducibly split into a scalar, 'J, which is the part of the vector parallel to n a , and a 
vector, $? a , lying in the sheet orthogonal to n a ; 

^ Q = ^n a + * Q , where * = ip a n a , and 4- Q = N ab ij) b = V> a , (9) 

where we use a bar over an index to denote projection with N ab . Similarly, any PSTF tensor, ip ab , can now be split 
into scalar, vector and tensor (which are PSTF with respect to n a ) parts: 

tpab = Tp(ab) = * (n a n b - \N ab ) + 2^ {a n b) + # a6 , (10) 

where 

<J ee n a n b i; ab = -N ab i; ab , 
* a ee iy a 6 n c ^6c = * 8) 

* Qb ee ^{afc} = (N (a c N b) d - ±N ab N cd ) ^ cd . (11) 

We use curly brackets to denote the PSTF with respect to n a part of a tensor. Note that for 2nd-rank tensors in 
the 1+1+2 formalism 'PSTF' is precisely equivalent to 'transverse-traceless'. 3 Note also that h{ ab } = 0, N/ ab \ = 

- n (a n b) = A ofc - \h ab . 



3 Our use of the term 'transverse' only refers to the fact that the tensor is orthogonal to n a \ this does not imply it is divergence free, 
which is an additional property of tensors in the commonly used TT gauge of the plane wave approximation. 
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We also define the alternating Levi-Civita 2-tensor 

e a b = £ab c n c = u d i] dabc n c , (12) 

so that e ah n h = = e iab)l e abc = n a e bc - n b e ac + n c e ab , e ab e cd = N a c N b d - N a d N b c , e a c e bc = N ab , and e ab e ab = 2. 
Note that for a 2-vector ty a , e ab may be used to form a vector orthogonal to , J but of the same length. 

With these definitions, then, we may split any object into scalars, 2-vectors in the sheet, and transverse-traceless 
2-tensors, also defined in the sheet. These three types of objects are the only objects which appear, after a complete 
splitting. Hereafter, we will assume such a split has been made, and 'vector' will generally refer to a vector projected 
orthogonal to u a and n a , and 'tensor' will generally mean transverse-traceless tensor, defined by Eq. {HJ. 

There are two new derivatives of interest now, which n a defines, for any object ip...'": 

i> a ... b c - d = n e D e W. b c - d , (13) 
6 e j, a ... b e " d = NjNj ■■■N b a N h c ---N i d T> j i>f... g h - i . (14) 

The hat-derivative is the derivative along the vector field n a in the surfaces orthogonal to u a . This definition represents 
a conceptual divergence from the tetrad approach, in which the basis vectors appear on an equal footing [i.e., with 
V a rather than D a in Eq. I|13f) ]. As a result, the congruence u a retains the primary importance it has in the 1+3 
covariant approach. (We choose to think of A = u a n b V a u b = —u a u h X? a n b as the radial component of the acceleration 
of u a , rather than the time component of fi a .) The ^-derivative, defined by Eq. (|14|l is a projected derivative on the 
sheet, with projection on every free index. 

With these definitions we may now decompose the covariant derivative of n a orthogonal to u a : 

B a n b = n a a b + \<j)N ab + ^e ab + Cab, (15) 

where 

a a = n c D c n a = n a , (16) 

<t> = S a n a , (17) 

£ = \e ab 5 a n bl (18) 

Cab = 5 {a n b} . (19) 

We may interpret these as follows: travelling along n a , </> represents the sheet expansion, Cab is the shear of n a 
(distortion of the sheet), and a a its acceleration, while £ represents a 'twisting' of the sheet - the rotation of n a [l7j . 
The other derivative of n a is its change along u a , 

n a = Au a + a a where a a = n a and A — n a u a . (20) 

The new variables a a , 4>, £, C, ab and a a are fundamental objects in the spacetime, and their dynamics gives us 
information about the spacetime geometry. They are treated on the same footing as the kinematical variables of u a 
in the 1+3 approach (which also appear here). 

Note that for a scalar, we have D a \l/ = $n a + <5 a \l/, while for any vector ^ a orthogonal to n a and u a (i.e., = ^ a ), 
we may decompose the different parts of its spatial derivative: 

D a * fc = -n a n b V c a c + n a i> b - n b [±0* a + (£e ac + Cac) * c ] + 5 a <$ b . (21) 

Similarly, for a tensor ^ a b- ^ab = ^{ab}, we have 

D a *6c = -2n a n {b ^! c)d a d + n a i> bc - 2n (b |^ c)o + ^ c) d (& ad + Cad) + 5 a ^ bc . (22) 

Note that N ab = 2u {a u b) - 2n {a n b ), N ab = -2n (a a b) , S c N ab = 0; while e ab = -2u [a e b]c A c + 2n [a £ b]c a c , i ab = 
2n^ a e b ] c a c , S c e ab = 0. 

We take n a to be arbitrary at this point, and then split the usual 1+3 kinematical and Weyl quantities into the 
irreducible set {6, A, SI, E, S, H, A a , S a , £ a , H a , S a6 , E abl H a b} using @ and ^ as follows: 

ii a = An a +A a , (23) 

uj a = nn a + fl a , (24) 

a ab = S (n a n b ~ \N ab ) + 2E( a n fc) + E ah , (25) 

E ab = £ (n a n b - \N ab ) + 2£ {a n b) + £ abl (26) 

H ab = H(n a n b -±N ab )+2H {a n b) +H ab . (27) 

Having described the splitting of the 1+3 variables to obtain their 1+1+2 parts, and the introduction of the new 
1+1+2 variables corresponding to the irreducible parts of V a rt , it only remains to apply this splitting procedure to 
the 1+3 equations themselves. We give these equations in section ITTT1 linearised about a Schwarzschild background. 







C. The Ricci identities 



Once the vector n a has been introduced it is possible, and necessary, to augment the 1+3 equations with the Ricci 
identities for n a ; without these we do not have enough equations to determine the new 1+1+2 variables. The Ricci 
identities for n° are 

J? a6c = 2V [o V 6] n c - Rabc d n d = 0, (28) 

where R a bcd is the Riemann curvature tensor. This 3-index tensor may be covariantly split using the two vector fields 
u a and n a , and gives dynamical equations for the covariant parts of the derivative of n a (namely a a , a ai <f>, £ and Cab) 
in the form of evolution equations, involving dot-derivatives of these variables, and propagation equations, involving 
hat-derivatives. In order to facilitate the calculation of these Ricci identities, which appear in the following section, 
we give here the expression for the full covariant derivative of n a in terms of the relevant 1+1+2 variables: 

V Q n b = -Au a u b - u a a b + (£ + \Q) n a u b + [E - e ac VL c ] u b + n a a b + \(f>N ab + £e ab + (at, (29) 

which may be inserted into Eq. (|28|l . 



D. Commutation relations 

In general the three derivatives we now have defined, ' ' ', and c S a ' do not commute. Instead, when acting on a 
scalar i/>, they satisfy: 

i-j> = -M+(^ + ^)^+{^a+e ab n b -a a )S a i;, (30) 
S a1 jj - N a b {8 b rl>y = -Aai> + (a a + E a - e ab Q b ) $ + (§0 - |E) 5 a ip + (Sah + fie„ 6 ) 5 b ijj, (31) 
5 a ^-N b (5^) = -2e ab Q. b iP + aJ+±4>6 a i> + (( ab +& ab )5 b iP, (32) 
6 a 5 b ip - 5 b 5 a il> = 2e a6 - £V>) + 2a[ a 5 b ]ip. (33) 

These last two equations are the decomposition of the 1+3 commutation relation 

cur]D V = 2^u; . (34) 

These relations are considerably more complicated for vectors and tensors. 

From Eq. (|33|l . we see that our sheet will be a genuine 2-surface in the spacetime (and, in particular, that the 
derivative 5 a will be a true covariant derivative on this surface) if and only if £ = = a a = 0. (Recall that the 1+3 
spatial metric h ab corresponds to a genuine 3-surface when w a = 0.) Otherwise, the sheet is really just a collection 
of tangent planes. In addition, the two vectors u a and n a are 2-surface forming if and only if the commutator [u, n] 
in 1)30(1 has no component in the sheet: that is, when Greenberg's vector, E° + e ab Q b — a a , vanishes [lij - see Eq. 1(30(1 . 



III. PERTURBATIONS AROUND A SCHWARZSCHILD BLACK HOLE 



For an exact Schwarzschild black hole it turns out the the only non-zero 1+1+2 variables are the scalars {^4, £ , 4>} 
(and their derivatives {A, £,(/>}), a covariant characterisation of the Schwarzschild solution. (We saw in section III Al 
that A and £ are the only non-zero parts of the 1+3 variables, and it is clear by considering Gauss' theorem that 0, 
the divergence of n a , must also be non-zero. We will consider the background solution in more detail in section Till Bh 
Because the background solution involves only scalars, under any perturbation all vectors and tensors are first- order, 
which greatly simplifies things, as we discussed in section III Al 

The usual 1+3 evolution and constraint equations may be further split with the vector n a , into a set of evolution 
(along u a ) and propagation (along n a ) equations. Together with the Ricci identities for n a , we find a complete set of 
propagation, evolution and constraint equations - the constraints being those equations with no hat- or dot-derivatives 
in them. We will give the complete nonlinear equations elsewhere, as they are large and unpleasant. Here, however, we 
will give the vacuum equations linearised around a Schwarzschild black hole background. Our linearisation procedure 
is straightforward: as in the 1+3 approach we neglect all products of first-order quantities; first-order quantities being 
those which vanish in the background: 



{9, n, S, £, H, a a , a a ,A a , fl a , E a , £ a , H a , Cab, Z ab , £ ab , H ab } = O(e) 



(35) 



(along with their derivatives, and dot- and (^-derivatives of {A, £,(/>}), where e is a 'smallness' parameter, which 
measures departures from an exact black hole. So, for example, one could define e ~ £l/4>, or e ~ \J £ ab £ ab /£, and 
so on. From now on all equations are linearised about a Schwarzschild black hole, and equations of the form A = B 
generally mean A — B = C(e 2 ) (in keeping with usual practice, we will not distinguish between this and real equality). 

When studying cosmological perturbations using the 1+3 approach, the evolution equations are of prime impor- 
tance, since time is the only remaining essential parameter: the goal is to find the evolution of seed perturbations 
corresponding to the various spatial harmonics. In contrast, for the black hole perturbations analysed here the time 
invariance and spherical symmetry of the background mean that radius is the interesting parameter, and so the 
propagation equations are the key: we want to find the variation with radius of the various (time and spherical) 
harmonic components. Thus, we present the propagation (hat) equations first, and relegate the evolution equations 
to a secondary position. This will be helpful when we come to solve the equations in Section Ivl 



Propagation 4 : 



-¥ 2 



Evolution 5 : 



-£ + S a a a , (36) 

i = -fc + \e ab 5 a a b , (37) 

a- a - a- a = - (y + A) a a + (±0 - A) (S Q + e ab fl b ) - e ab H b , (38) 

|(9-E = §0£ + <5 Q £ a +e ab d a n b 1 (39) 

A- 6 = -A{<t> + A) -5 a A a , (40) 

d = {A-<p)n-5 a n a , (41) 

£ = -<5 a £ a , (42) 

n = -^n-5 a n a -sm, (43) 

E a -e ab Cl b = -l^ a + (2A+±cb)e ab n b + lS a 6+±8 a Z~e ab 5 b n~5 b Z ab , (44) 

A- a -2±- a = 2£ a — Aa a — S a A — 2 (A — \<f) A a , (45) 

4 = + \5 a £ - \£ a a - S b £ ab , (46) 

H- a = -^H a + ±S a H-S b H ab + l£(n a -£ ab Z b ), (47) 

E{ab} = -\^ah + 5{ a ^ b } - £c{a& C Q b } - £c{a7~(. h } , (48) 

C{afc} = -4>(ab - £ab + S{ a a b }, (49) 

£ {ab} - e c{a H b} c = -l£Z ab + e c{a (-6 c H b} + (^ + 2A)H b} c ), (50) 

n {ab} + e c{ J b} c = e c{a (d c £ b} + |£C 6} C - (\4> + 2A) £ b} c ) . (51) 

4> = {\<t>-A) (E- \9) +6 a a a , (52) 

£ = -(\(t>- A)VL + \H + \e ab 5 a a b , (53) 

tl = £A+h £ abS a A b , (54) 

\9-t = £ + <t>A + $aA a , (55) 

£ = (^-e)£ + e ab S a H b , (56) 



4 These equations are derived as follows: Eq. 1361 from n a N bc R a i, c ; Eq. 1371 from n a e bc R a b c \ Eq. 1381 from u a n b R. a i,c = n a u b R a bc\ 
Eqs. 1391 and 1441 from the shear divergence equation, (C\) a \ Eq. 1401 from the Raychaudhuri equation; Eq. 1411 from the rotation 
divergence equation, C2; Eqs. 1421 and 1461 from the electric Weyl divergence equation, (Cj)"; Eqs. 1431 and 1471 from the magnetic 
Weyl divergence equation, (C5) 11 ; Eq. 1451 from the shear evolution equation; Eq. 1481 from the iT a j,-equation, (C'i)^ a i,y, Eq. 1491 from 
n<1 Ra{bc} \ Eq. I5UI from the electric Weyl evolution equation; and Eq. 1511 from the magnetic Weyl evolution equation. 

5 These are derived as follows: Eq. 1521 from u a N bc R a bc\ Eq. 1531 from u a e bc R a b c \ Eq. 1541 from the rotation evolution equation; Eq. 1551 
from u a n b u c R a b c = —n a u b u c R a b c and the Raychaudhuri equation; Eqs. 1561 and 1591 from the electric Weyl evolution equation [ 1591 
also uses l47l l ; Eqs. 03 and ^ from the magnetic Weyl evolution equation f i60l also uses 1461 1: Eq. 1581 from the rotation and shear 
evolution equations; Eq. 1611 from the shear evolution equation; and Eq. 1621 from ti c _R c { ai) j. 



Constraint 6 : 



H = -3£^-e ab S a £ b , (57) 

t- a -e ah tt = -£ a + 5 a A+(A-\4>)A a , (58) 

4 = -f £oa + |e ab {5 b U +{4>- 2A) H b ) - s c{d 5 d H a f, (59) 

Ha = -§£ £ „ b .A b - f £ o6 (5 6 £ + (<(>- 2A) £ b ) +e c{d S d £ a} c , (60) 

£{ah} — -^Cafc - £ab + 5{ a A b }, (61) 

C{afc} = ~ {\4> ~ A)Y, ab - e c { a U b f + 8{ a a b }. (62) 



(S„E - §<5 a £? + 2e ab 5 b n + 26 b E ab = —<j>Yl a + <ps ab Q b — 2s ab H b , (63) 
6 a n a +e ab 6 a Z b = (2A - 4>) + W, (64) 
§<W - e«6* 6 £ - <* 6 C«& = (65) 

These equations must of course be consistent with one another: the constraints must evolve and propagate consis- 
tently, and each first-order variable must satisfy the commutation relation 

C[$] = # - I - AV = 0, (66) 

while the background scalars must satisfy the commutation relation Q3U[I. It turns out after an arduous calculation 
that all the equations are consistent. 



A. Frame choice 



In this work we are presenting a gauge invariant, covariant approach to perturbations of spherically symmetric 
spacetimes, based on the introduction of a partial frame, that is, of two basis vectors. In GR there are two 'gauge' 
freedoms: the choice of coordinates and the freedom to choose a frame basis in the tangent space at each point. These 
have their direct analogues in perturbation theory, where we imagine that the true spacetime we are studying is 'close 
to' some given, idealised background spacetime. Although not formally the same thing, choosing a coordinate system 
in the true spacetime is in practice equivalent to fixing the mapping between the true and background spacetimes 
allowing the direct comparison of scalar, vector and tensor objects in the two spacetimes at corresponding points 
(see [13). In metric- based (non-covariant) approaches to the perturbation problem the first step is to find a nice 
coordinate system in the true spacetime, corresponding to that in the background, and to write equations for the 
derivatives of the perturbations of scalars, vectors and tensors with respect to these coordinates. 

In covariant (partial-)frame approaches, on the other hand, one tries to avoid explicit reference to the background, 
using it merely to determine which quantities are zeroth order (i.e., which do not vanish in the background). Given 
the frame vectors, a set of covariant (that is, coordinate invariant) equations describing the true spacetime are written 
down. (Coordinate-)gauge invariance will hold, according to the Stewart- Walker lemma if we can find a complete 
set of variables all of which vanish in the background. However, since the covariant variables are the projections of 
tensors with respect to the frame vectors and the projected parts of the of frame-vector covariant derivatives, the 
equations and their solutions will, in general, depend on the particular choice of frame vectors. However, since the 
true spacetime lacks the symmetry of the background, there is, in general, no unique covariant definition of the frame 
vectors, and one can always perform a first-order rotation of these. This freedom can easily be seen in the preceding 
set of equations: there are no evolution equations for A, A a , and a a , while there is no propagation equation for a a . 
Indeed, this is true in any spacetime, as one can choose the frame vectors at any point freely, the motion of which 
must be put into the equations by hand (GR can't predict this!). 

In what follows we will reserve the term 'gauge invariant' to refer to the invariance of the equations under the 
mapping between the true and background spacetimes (in the sense of the Stewart- Walker lemma), and will use 
'frame invariant' to describe invariance under the choice of frame vectors. Our 1+1+2 formalism, in common with 
the 1+3 and Newman-Penrose approaches, is gauge invariant, but not frame invariant. 



These are derived as follows: Eq. 1631 from (Cz) ab n b and {C\)a, or ri a u c R a i c ; Eq. 1641 from e ab u c R a f >c ; and Eq. 1651 from N bc Rabc- 
Note that the equation formed from (Cs) a i,n a n h is equivalent to Eqs 1411 and 1641 . 
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The equations as they are presented above are completely general, involving no particular choice of either frame 
vector. Given a timelike congruence u a , obvious choices are possible: for example, we could take n a to be parallel 
to the acceleration u a , leading to the frame condition A a = 0, a physically plausible choice for observers hovering 
above the black hole as one could always dexterously align ones rocket to make this so (objects dropped onto the 
floor of the rocket would fall directly 'down'); or we could demand that n a be parallel to the 'radial' eigenvector of 
the electric Weyl tensor E a b, which leads to the frame choice £ a = 0. Indeed, we can set any 2- vector to zero by these 
considerations. We shall not impose a frame condition at the moment, however, as there is no need. 



B. Background Solutions 

In the background we have, by setting all vectors, tensors and time derivatives to zero and retaining only the 
zeroth-order scalars: 

$ = -\4> 2 -£, (67) 
£ = ~l4>£; (68) 

together with 

£ + A(f> = 0, (69) 

which is Eq. Q55[l. Note that these three equations are sufficient: Eq. (|40|l is satisfied. 

If we associate our hat-derivative with an affine parameter p, i.e., * = d/ dp, then we may solve these equations. 
The parametric solutions to our background equations are, in terms of either a parameter x or r, 

r sech 6 x = -, (70) 



(2m) 2 r 3 



1 , 2 / 2m 

6 = — sinhrrscch x—-\ I , (71) 

to r V r 

1 o m ( 2m\ 1//2 

A = - — cosech x sech 1 x = — I 1 ; (72) 

4to r z \ r I 



where 



and the usual Schwarzschild coordinate 



2m {x + sinh x cosh x} , (73) 



= 2m cosh 2 x. (74) 



These form a one-parameter family of solutions, parameterised by the constant to, which is just the Schwarzschild 
mass. The Schwarzschild solution is given for 2m < r < oo for < x < oo; the interior solution may be found by this 
approach, but we will not require it here. 

We show a plot of <\> and A with r in Fig. ^ This shows how the expansion of n a starts from zero at the horizon, 
is largest at the photon sphere, before dropping to zero again as r — > oo. 

The solutions given by Q72JI represent the general solution to the system of equations l|67|l. I|68(l and H69J) . which 
are the covariant 1+1+2 equations under the conditions of the spacetime being static, spherically symmetric and 
vacuum; hence, we may covariantly characterise the Schwarzschild solution by {A, £,4>} ^ 0, with all other covariant 
quantities being zero. 7 



C. Gauge invariant variables 

We have developed a set of covariant equations describing a perturbed black hole, with all variables defined with 
respect to a family of observers u a with a preferred radial vector field n a , which these observers can choose; first- 
order variables have a clear physical or geometrical meaning. However, not all variables appearing in the equations 



7 We would like to thank an anonymous referee for drawing this to our attention. 
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FIG. 1: A plot of (j> and A with r, showing the maxima of <j> at the photon sphere, r = 3m. In contrast, A falls from oo at the 
horizon. 

are gauge invariant, because they are not first order. Recall the Stewart- Walker lemma 0|j which states that if a 
variable vanishes in the background then it is gauge invariant in the perturbed spacetime. 

Equations (f4*n)l . (|55|l . and arc not gauge invariant because zeroth-order terms appear in these equations, 
i.e., isolated terms involving {A,£,4>} not multiplied by a first-order quantity. We therefore introduce the set of 
gauge-invariant variables 

X a = 5 a £, (75) 
Y a = <5a0, (76) 
Z a = 5 a A, (77) 

which do vanish in the background ({A, £, <j>} do not vary over a sphere in the background), and so are gauge invariant, 
by the Stewart- Walker lemma 0] . It is important to notice that we lose no true degrees of freedom in the solutions 
to the equations by introducing these variables, since we only eliminate possible spherically symmetric perturbations 
(for which X ai Y a and Z a are automatically zero), but we know from Birkhoff's theorem that all spherically symmetric 
vacuum spacetimes are Schwarzschild, and therefore any such nontrivial solution must be purely a result of the freedom 
to chose the frame vectors in the Schwarzchild background. 

We may directly obtain the following evolution and propagation equations for these new gauge- invariant variables: 



x a = ys(a a + i: a -£ ab a b ) + ^£(s a i:-y a e)+£ bc s a s b n c , (78) 

Y a = (i<^ 2 +£)(a a + Y, a -e ab n b ) + (\4>-A)(8j:~l5 a e)+5 a 5 c a c , (79) 

X a = -2(bX a - \£Y a + %cj>£a a - 5 a S b £\ (80) 

% = -X a -'^Y a + (^ 2 +£)a a + 5 a S b a b , (81) 

Z a = -(l0 + 2A)Z a -AY a +A(^ + A)a a + 5j-5 a 5 b A b . (82) 

Note that there is not an equation for Z a because there is no equation for A. 

We also find the following constraints by applying the commutator to our new variables: 

e ab 8 a X b = 30£C, (83) 

e ab S a Y b = (^ 2 + 2£)£, (84) 

e ab 6 a Z b = 2A(4> + A)£. (85) 



Equations l|78|l and l|79J) replace our equations for H56(l and l|52|l respectively. Similarly Eq. I|8()(l replaces Eq. H42J) . 
Eq. (|81|l replaces , and (|82|l replaces (|4(J|I . The constraint l|85|) allows us to derive a pseudo-evolution equation for 
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Z a (but we shall not require it). We may also replace our £ equation l|55|l with 

6 a ± - §<y = -X a - AY a - ^Z a - 8 a S b A b . (86) 

All equations are now gauge invariant. This means that when we have chosen our frame u a ,n a in a unique way, 
then all quantities appearing in the equations are uniquely defined covariant and gauge-invariant quantities, with a 
physical or geometrical meaning. 



D. Commutation relations for first-order variables 



Having discussed which variables are zeroth order and which are first order, it will be useful to present the following 
commutation relations for the derivatives of first-order scalars, vectors and tensors. For any scalar, vector or tensor, 

V-4> = AV, (87) 
5 a 4> - (S a ^y = 0, (88) 
S a 4f-{S^) = U5 a V; (89) 



while for a scalar 



and for a vector 



and a tensor 



5 [a S b] * = 0, (90) 
6[a6 b ]* c = (\<t> 2 -£)N c[a * b] , (91) 
5 [a S b] f c d = \ (\4> 2 - £) [N c[a V b]d + N d[a V b]c ] . (92) 



E. Spherical harmonics 

As the equations stand we can't find solutions because of the appearance of angular derivatives : S a \ An appropriate 
choice of basis functions will allow us to write all first-order variables as an infinite sum over these basis functions, and 
allow us to replace angular derivatives by a harmonic coefficient. Clearly the spherical symmetry of the background 
begs us to use spherical harmonics as our basis functions, so here we will develop these appropriately for our formalism. 
We define our harmonics by analogy with the FLRW case [|, and we refer to |2(j for details of other approaches. 

Note that all functions and relations below are defined in the background only; we only expand first-order variables, 
so zeroth-order equations are sufficient. 

We introduce spherical harmonic functions Q = Q^> m \ with m — —£, ■■-,£, defined on the background, such that 

6 2 Q = -£{£+l)r- 2 Q, Q = = Q. (93) 

The function r is covariantly defined by 

= 2-, f = = S a r. (94) 
r 

This factor is included in our definition (|93|l so that the equation propagates; it is trivial to show that it evolves also. 
We have defined r so far up to an arbitrary constant, which reflects our freedom in choosing a particular normalisation 
of the spherical harmonic functions; we will find it most useful for our purposes to fix this freedom by covariantly 
defining 

r=(^-£)- 1/2 , (95) 

i.e., we identify r defined here with the parameter defined by Eq. H74|) . We can now expand any first order scalar \t 
in terms of these functions as 

oo m—£ 

* = J2 Yl * { s ,m) Q ie ' m) = *sQ, (96) 

1=0 m =-l 
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where the sum over £ and m is implicit in the last equality. The S subscript reminds us that ^ is a scalar, and that 
a spherical harmonic expansion has been made. Due to the spherical symmetry of the background, m never appears 
in any equation; we will just ignore it from now on. 

We also need to expand vectors and tensors in spherical harmonics. We therefore define the even (electric) parity 
vector spherical harmonics for £ > 1 as 

Q [ a ] =rS a Q {l) Qa = = Q a , S 2 Q a = (1 - £ (I + 1)) r- 2 Q a ; (97) 

where the {£) superscript is implicit, and we define odd (magnetic) parity vector spherical harmonics as 

Q { a e) =re ab S b Q^ =► Q a = = Q a , S 2 Q a = (1 - £ (£ + 1)) r - 2 Q a . (98) 

Note that Q a = e a bQ b & Qa = —£abQ b , so that e a b is a parity operator. The crucial difference between these two 
types of vector spherical harmonics is that Q a is solenoidal, so 

S a Q a = 0, while S a Q a = -£(£+ l)r _1 Q. (99) 

Note also that 

£ ab S a Q b = 0, and e ab 8 a Q b =l(£+l)r~ 1 Q. (100) 

The harmonics are orthogonal: Q a Q a — (for each £), which implies that any first-order vector \l/ a can now be 
written 

oo 

*a=E ^VQ^ + *W = *VQ Q + * V g Q - (101) 

1=1 

Again, we implicitly assume a sum over £ in the last equality, and the V reminds us that ^ a is a vector expanded in 
spherical harmonics. 

Similarly we define even and odd tensor spherical harmonics for £ > 2 as 

Qab = r 2 5 {a S b} Q, => Q ab = = Qab, 6 2 Q ab = [4> 2 - ZE - £ {£ + 1) r- 2 ] Q ab ; (102) 
Q ab = r 2 e c{a S c 6 b} Q, => Q ab = = Q ab , S 2 Q ab = [4> 2 - 3£ - £ {£ + 1) r" 2 ] Q ab , (103) 

which are orthogonal: Q a bQ ab = 0, and are parity inversions of one another: Q ab — —e c {aQb} ^ Qab — £ c {aQbf- 
Any first-order tensor may be expanded 

oo 

*ab = ^Qab + *T ] Qab = *tQ«4 + *jQab- (104) 

1=1 

We will not write all the equations expanded in spherical harmonics; instead we list here all the replacements which 
must be made for scalars, vectors and tensors. Note that sums over £ and m are implicit in these equations below. 
For brevity, we will sometimes use the aliases 

L = £{£+!), I = {£- 1) (£ + 2) = L - 2. (105) 



SCALAR VECTOR TENSOR 



# = 








*a& = 


+ ^jQab 


+ ^jQab 


Sa^ = 


r-^sQa 


£ab^ b 


= -*V<9a + *vQa 




-^jQab 


+ ^jQab 


£ab8 b ^ = 






= -£(£+l)r- 1 V v Q 


S b ^ab = 


\lr-± (- 


^TQa + ^jQa) 






£abS a ^ b 


= +£(£+l)r- 1 ^ y Q 


£c {d 5 d *: } = 


fr- 1 (+ 


#T<9a + ^jQa) 






S{a^b} 


= r" 1 (* V Qab - *vQafc) 












£ c {J c ^b} 


= r- 1 (* V Qa6 + *vQ Q b) 









Each vector and tensor equation produces two harmonics equations for each £, one of odd parity and one of even 
parity, due to the orthogonality of the vector and tensor harmonics. 

We can use these harmonic relations to derive various properties of the 6 a derivative. For example, an important 
relation we use in deriving the Regge- Wheeler equation in Sec. HVl is 



25 {a 5 C ^ b }c ~ S 2 ^ ab =(£- *ab- 



(107) 
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F. Odd and even parity perturbations 

After decomposing the equations into their harmonic components, we find that they split into two independent 
subsets, which we refer to as even and odd (parity) perturbations, but are also known as polar and axial perturbations. 
This splitting is analogous to the well known splitting of scalar, vector and tensor modes in the cosmological situation. 
The two independent sets of equations utilise the variables 

Odd parity: V = ({£ T , H T , S T , Ct}, {£v, Hv, Ev, fiv, Ay, ay, ay, Xy, Yy, Zy}, {H s , fi S , &})(108) 
Even parity: V E = ({£t, Ht, St, Ct}, {£v, Hv, S v , fiv, Ay, ay, ay, Xy, Yy, Zy}, {S s , S }) , (109) 

where we have defined the vectors Vo, Ve for later convenience. The 'parity switching' which occurs between the sets 
of variables (e.g., Hj appears in the odd parity system) may be seen in the covariant tensor equations: H a b always 
appears alongside a l £ a b factor relative to other variables such as £ a b\ similarly for fi a ; it's always seen as s a b fi 6 , 
relative to £ a , say. 

Hereafter, we will assume all equations have been decomposed into their spherical harmonic components, unless 
stated otherwise, and when we refer to specific equations that are given in tensor form (such as the evolution equations 
above), we will generally assume too that this has been decomposed implicitly. 



G. Time harmonics 



Because the background is static, we can, if we wish, decompose time derivatives of first order quantities into their 
Fourier components. This is simply assuming an e lulT time dependence for the first order variables, but we shall make 
it a bit more precise. Define the time harmonic functions in the background by 

fM=iwT( u \ f ("> = = S a T^; oj = = 5 a u. (110) 

In the background, this must satisfy 

T + ^T = 0, (111) 

which implies 

u> = -Au. (112) 
We may integrate this in the background, in terms of the parameter x, or r, to give 

/ 2mV V2 2a . „, 

w = crcothx — all = — , (113) 

\ r J 4> r 

where a is a constant. Then any first order variable ^ (which will usually be an even or odd parity variable discussed 
above, but it could be any first order scalar, vector or tensor) as 

v]/ = vJ/MtM — xj/MyM (114) 

where the summation is understood in the last equality (which may be an integral, and depends on the types of 
differential equations and their boundary conditions occurring in the solutions). We can simply replace a dot by iu) 
in the equations, as no confusion should occur. 



IV. THE REGGE-WHEELER TENSOR AND THE REGGE- WHEELER AND ZERILLI EQUATIONS 

We have presented the full, covariant, gauge-invariant linearised equations for the propagation of gravitational 
radiation in a perturbed Schwarzschild black hole spacetime, and we have discussed the introduction of spherical 
harmonics, enabling us to replace ^-derivatives with spherical harmonic indices, as well as remove the tensorial nature 
of the equations. We could now introduce a perturbed metric and calculate all variables in terms of the metric 
functions, to show how the standard Regge- Wheeler and Zcrilli equations of black hole perturbation theory may be 
related to this approach, linking all of our variables to these Regge- Wheeler and Zerilli functions. We need not bother 
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however: it is possible to find generalisations of these functions and correlations directly, an important test of our 
theory. 

We show in this section that it is possible to find a gauge- and frame- invariant, transverse-traceless tensor that 
satisfies a closed, covariant, gauge-invariant wave equation. We call this tensor the Regge- Wheeler tensor. Moreover, 
we demonstrate that, once decomposed into spherical harmonics, and with the appropriate radial coordinate, this 
equation is the Regge- Wheeler equation for both parities, thus unifying both parities in one TT tensor. We present 
similar results for the Zerilli equation that normally describes even parity solutions, but defer the derivation of the 
Zerilli tensor until later in Sec. IV CI 



A. Regge- Wheeler 

We know that gravitational waves propagate in the perturbed Schwarzschild spacetime, and therefore that the TT 
tensors must satisfy (covariant) wave equations of some sort [see Eqs. <|152|) and (|153[) in section lVB| . as our solutions 
must recover the plane wave case when m = 0. When we investigate these wave equations for £ ab and Ti a b, say, we 
find that they are not closed: that is, they contain forcing terms from other 1+1+2 tensors, a feature not present for 
plane waves. This makes their interpretation and solution nontrivial. In fact, when we look at the second-order wave 
equations obtained from the covariant equations for any of the 1+1+2 tensors we find the same story. Can we find 
some combination of the basic tensors that satisfies a closed, covariant wave equation? Indeed we can, and we outline 
its derivation here. The formulae are rather formidable, and so we omit a detailed derivation. 

In addition to the obvious 1+1+2 transverse-traceless tensors (£ a b, "Hah, £a& and Cab) it is possible to construct 
many TT tensors from (5-derivatives of vectors and scalars, such as 5{ a X b y, 5{ a a b y, or even S[ a 5 b yfl, for example. 
Using the dot and hat equations, along with the commutators, we can obtain the wave equations satisfied by any of 

these tensors, by calculating the wave operator ^ { a b} ~ *{a&} f° r that tensor ^ ab . 

If we calculate the wave operators for Cab and 8{ a X b y we notice that they contain similar terms in a a , a a and (-. 

In fact, if in ({ a b} — C{ a b} we substitute for a a from Eq. JHOJl, a„ from Eq. (|7H|l . and for £ from Eq. lf%3|) we find, 
amazingly, that all variables other than 5{ a X b y and £ ab cancel, leaving the closed wave equation (|ll(j[) below. 
So, we define the dimcnsionless, gauge-invariant, frame-invariant, transverse-traceless tensor 

W ab = ^r 2 ( ab -±r 2 £- 1 5 {a X b} , (115) 

which obeys the rather nice wave equation 

W {ab} - W {ab} - AW {ab} + tfW ab - S 2 W ab = 0, (116) 

where S 2 = S a S a is the covariant laplacian on the sheets (approximate 2-spheres, in this case). 
We can expand Eq. (|TTE|) into SH: let W = {W T , W T }, then Eq. i|TiT)Jl becomes 

W-W- AW + ( - ^ t + 3g \ w = °- 



Note that both the even and odd parity parts of W ab satisfy the same wave equation (I117J) . 

It turns out that Eq. (|117fl is actually the Regge- Wheeler equation |8| in appropriate coordinates, which we now 
show. Converting to the parameter r,p^>r, and then to the 'tortoise' coordinate of Regge and Wheeler 



and then letting 

we find that (|117|l becomes: 

where 



r + 2mln(— - lV (118) 
\2m J 



^ = <pRW = W, (119) 



-^ + CJ 2 jr = lV (120) 



V = V RW = {r . 2m) (l(l+l)r-6m), (121) 



which is the Regge- Wheeler potential. We will thus refer to W ab as the Regge- Wheeler tensor. 
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B. Zerilli 



We have shown that W a b satisfies the Regge- Wheeler equation (|117|) regardless of parity. However, it is well known 
that even parity perturbations are governed by the Zerilli equation. 
For even parity perturbations the variable 

Z = ^{l-ZEr 2 )^ 1 (3rZ T + 2£- 1 Hv) (122) 
[3r0S T - 2A~ 1 H\,] (123) 



6 
3 C 3 



is one of two fundamental variables, and it can be shown to satisfy the Zerilli equation. The wave equation equation 
for Z is 



\r- 2 { ±c 3 + 32C3- 2 (L + 1) I 2 } Z = 0. (124) 



We have used the abbreviations: 



c,=4(L + 1)-jW ! (125) 



where we keep the freedom in j for later use. 

Making the change to the tortoise coordinate, inserting time harmonic functions, and changing to the variable 

i/> = = Z (126) 

we find that ij 124(1 is in fact the Zerilli equation, Eq. 1(120(1 . with 

2 "' \{L- if [r 3 (L + 7) + 24mr 2 ] + 36m 2 r (L - 1) + 72m 3 ] . (127) 



V = Vz = 

r 4 {(L- l)r + 6to} 2 - 

For the even parity perturbations, then, there are two variables which obey wave equations: the Zerilli variable Z, 
and the Regge- Wheeler variable Wj. They are in fact related by 

luJ r 2 2 [4 -8l(c 3 + L + 4)] + 32LI (L + 1) 
Z= ^A Wj + (18) 



as may be found by considering considering Wj, using Eq. 1(115(1 . Furthermore, it is possible to find an equation of 
the form 

Wj = (stuff) Z + (more stuff) Z, (129) 

by utilising Eqs. ((59(1 . ((50(1 and 1(61(1 . using the time harmonics throughout, which implies that Z may be written as 
a function of Wj and its derivatives. This gives 

. r 2 2 [ C 2 -8Z(c3 + £ + 4)]+32^(£ + l) w QuA^M 2 + L 2 Z% n _ 

Wj = w T sT^w — z - (130) 

Thus, the wave variables are not independent. Eqs. ((128(1 and 1(130(1 are a two dimensional first-order linear system 
of differential equations, which could replace the second-order Regge- Wheeler and Zerilli equations. In fact, it is 
instructive to rewrite Eqs. ((128(1 and ((130(1 in matrix form: 

(131) 

where the definitions of A, B and (3 are obvious from above. Notice that the matrix that couples Wj and Z is 
traceless, and that the special quasinormal modes discussed below correspond to A — 0. 

This form of writing the Regge- Wheeler and Zerilli equations, Eq. I(131|) . may be viewed as just a neat way to 
write two decoupled second order DE's as two coupled first-order ones; however, it also shows that the two decoupled 
second order DEs, which Eq. ((131(1 may replace, must be representations of the same physical situation (because 
Wj is a linear combination of Z and Z : and Z is a linear combination of Wj and Wj). We can use this DE to find an 
odd parity Zerilli variable, which satisfies the Zerilli equation ((124(1 . quite easily: simply demand that Eq. 1(130(1 hold 
for the odd parity perturbations too, and substitute for Wj from Eq. ((115(1 . and then from the propagation equations 
as appropriate. This then gives us Z as a complicated linear combination of £j, (j, Xy and Yy, seemingly unrelated 
to the even parity Zerilli variable, Eq. ((123(1 . This is an illusion, however, and it is possible to find a form for Z very 
similar to Eq. I(123|) . and consequently a Zerilli tensor, but we defer this until later, in Sec. IV Cl 
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C. Quasinormal modes 

The Regge- Wheeler and Zerilli equations have been studied in some detail over the years, and their solution is a 
complicated business [2l]. The relevant boundary conditions for the two equations are those that represent a GW 
perturbation which propagates outwards at infinity (r ~ r* — ► oo) and inwards to the horizon (r — > 2m, — * — oo) - 
that is, there are no GW propagating in from infinity or out of the horizon. The form of the Regge- Wheeler and 
Zerilli variables corresponding to this are 

ip ~ e lcrr * as r* — > — oo and tp ~ e _lcrr * as — > +00, (132) 

where ip — W or tp — Z; see, e.g., [2]]]. It turns out that the only solutions to Eq. (|120|l . with potentials 1(121(1 or (|127|) 
with boundary conditions i(132|) require discrete values of the frequency parameter a, with 5(cx) > 0; these are referred 
to as quasinormal frequencies, and the solutions constructed from them as quasinormal modes (QNMs). Because of 
the e lulT time dependence, these decay exponentially in time, which corresponds to energy radiated to infinity or the 
horizon as GW. This damping in time is important as if) grows exponentially as r — > 00. In particular, it also means 
that the spacetime is not flat at spacelike infinity, but it is flat at future null infinity (i.e., along a null ray). 
The factor in front of the Z in Eq. \ 130(1 is rather interesting: evaluate (in terms of r say) to find 

9uj 2 r 8 (j) 4 A 2 + L 2 l 2 = (12toct) 2 + (LI) 2 (133) 

which has roots at 

2mu = ±^£ {£-!){£ +!){£+ 2) (134) 
6 

with the '+' root corresponding to the frequency of the 'special (quasi-normal) mode' discussed in [2lt |22|. which is 
the only QNM with Ji(er) = 0. Since the frequency is purely imaginary the special QNM is not oscillatory in time but 
merely decays exponentially. 

The potentials for the Regge- Wheeler and Zerilli equations when converted into Schrodingcr form, Vbw and Vz, 
although appearing very different functionally, are actually very similar (see, e.g., 2 lj ) , becoming identical as £ — > 00, 
with peaks lying just beyond the photon sphere, r = 3m: as £ — > 00 the peaks approach the photon sphere (so that the 
solutions W and Z become identical for t — > 00). Thus these one dimensional wave equations represent gravitational 
waves scattering off the photon sphere, with the same reflection and transmission coefficients 0. 

It is important to note, however, that the Regge- Wheeler and Zerilli equations allow physical solutions with bound- 
ary conditions other than ((132(1 . and hence with ^s(a) < 0, but these solutions will represent GW incoming partially 
from infinity. These are not relevant for GW detection, but are of interest in their own right, as perturbations of 
bounded regions, say, may involve such waves. 



V. SOLVING THE EQUATIONS 



Before we discuss the solution of our system of equations, it is worthwhile giving an overview of the structure of 
the equations. We have three distinct types of equations: propagation, evolution, and constraint. The propagation 
equations are to be considered the key differential equations, while the evolution and constraints may be considered 
as auxiliary algebraic equations. This is because the structure of the background varies in the radial direction, so the 
hat-derivative cannot be expanded in harmonic functions, while the dot- and ^-derivatives are perturbation derivatives, 
as they do not occur in the background equations, and can thus be expanded in harmonics. We can then analyse and 
solve the large system of equations using matrix methods. This will provide interesting insight into the problem of 
black hole perturbations, since it allows us, at any radial position from the black hole, to treat the harmonic variables 
in (|108fl and (|109fl as 'coordinates' - i.e., as a particular choice of basis vectors - in an abstract, 33-dimensional 
vector space V33: the non-propagation equations then tell us that only fourteen degrees of freedom are present in 
the evolution and constraint equations (six in the odd parity, eight in the even), and these then allow us to work 
out how a subset of this this 14-dimensional solution subspace of V33 propagates radially (that is, they provide the 
differential system for the remaining degrees of freedom, but there are some undetermined variables), without us 
having to explicitly specify a basis for the 14~ dimensional vector subspace of V33. 

Assume a spherical harmonic decomposition of all the equations. Let V denote the 33-dimensional vector (element 
of V33) with the variables Vq = 1(11)8( 1 and Ve — i|lU9|) as elements; arrange the vector thus: 



V = (Odd variables | Even variables) = (Vq, Ve). 



(135) 
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Further, introduce the harmonic expansion in time, so that dot derivatives are everywhere replaced by iu)\ the evolution 
equations without a propagation derivative in them then become a set of 18 algebraic equations. The propagation 
equations are then of the form 

V 28 = PV, (136) 

where V28 denotes the vector consisting of the 28 elements of V which have a propagation equation. Note that there 
is no propagation equation for a a (indeed, a a only acts as a forcing term in the differential equations - it is completely 
undetermined) , and some variables do not have an individual propagation equation: there are propagation equations 
for £ a — e a b£l b and £ — but not for each of these variables separately (and these can't be generated, as the evolution 
equations and constraints contain the same degeneracy). P is the 28 x 33 propagation matrix, which contains iui 
terms subsuming evolution equations which have hat derivatives in them. 
The 18 remaining evolution equations take the form, in matrix notation 

EV = 0, (137) 

where E is a 18 x 33 matrix. The 9 constraints are 



CV = 0, (138) 

where C is a 9 x 33 matrix. Recall that the equations decouple into two sets of equations of opposing parity, which 
is reflected in the structure of the matrices E, P, C, which are all in block-diagonal form: 



Odd 


\ 


parity 





Even 


parity / 



thus dividing V33 into two divorced subspaces. (We will denote the odd parity upper block of the matrices by a 
subscript 0, and the even lower block by E.) This means that we can treat the odd and even subspaces separately. 
We will find that the odd equations reduce to a third-order system (i.e., a coupled three-dimensional first-order system 
of differential equations), and the even equations to a fourth-order system, both of which do not close: the odd has 
two undetermined variables (from the £ a — e a b^ b degeneracy, and lack of an equation determining a a ), while the even 
has three (from the £ a — e a b^ b and £ — |0 degeneracies, as well as the a a business), all of which reflect freedom in 
choosing our frame. 

In principle, we could adjoin the matrix C to E giving a total of 27 algebraic relations between the variables, 
represented by the 27 x 33 matrix 




(140) 



with 



FV = 0. (141) 

However, only one constraint equation, Eq. JSSJ, is not represented in the evolution equations (as there is no equation 
for Z a ). With this exception the constraints evolve and propagate consistently, so eight elements of C do not contain 
any more information than E and P. To put this more elegantly, the fact that the constraints evolve consistently 
implies that eight rows of C are just linear combinations of the rows of the matrix E, so that the rank of F is 19, 
just one more than that of E. Let Fc denote the 19 linearly independent rows of F. Formally, this means that 
the solution vector V in V33 must lie in the (33 — 19 = 14-dimensional) null space, A/14, of F^. Since all equations 
propagate consistently we can think of this solution space, and the propagation equations acting on it, in an abstract 
way. We see that the propagation equations are in a certain sense independent of the particular variables we choose 
to represent the solution. 

To actually obtain solutions, however, we must reintroduce 'coordinates' in A/14. To this end, we use F^ to eliminate 
19 variables, leaving just 14 (it being largely a matter of choice exactly which 14, provided they span A/14); let us 
denote them by v. The remaining variables can then all be expressed in terms of these fourteen 'coordinate variables'. 
We may therefore write 



V = Mv, 



(142) 
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where M is a 33 x 14 matrix of the form 

/ 



ODD 



V 



EVEN 



(143) 



/ 



Because of the frame freedom evident in the degenerate variables, we will split the vector v into two parts: v = 
(vo,vp), a 'determined' part containing 9 variables which have an individual propagation equation and a 'free' part 
containing 5 [say, (^v,av | ^v,av,#s)] which do not. Inserting Eq. 1(142(1 into the propagation equation, Eq. (|136fl . 
will then finish off the solution, resulting in the equation 

v D = Bv D +Av F , (144) 

where B is a 9 x 9 matrix and A is 9 x 5, giving the solution in the form of a nine-dimensional, first-order (non- 
autonomous) dynamical system for the vector ~v Dl with forcing terms from the five undetermined variables. The 
solution (114211 is then guaranteed to propagate using ^ 144(1 . 

From a naive argument based on counting degrees of freedom of our frame vectors we would expect to be able to 
eliminate five variables through a careful use of all frame freedom: the congruence u a can be changed by boosting 
to any new frame moving with some (first-order) three- velocity, giving three degrees of freedom; given a u a , we may 
further make any first-order change in n a that preserves u a n a = 0, giving two degrees of freedom. Indeed, there are 
five unknowns in Eq (1 144(1 , corresponding precisely to these five degrees of freedom, which implies that five equations 
(rows) of Eq. 1(144(1 do not represent true dynamical degrees of freedom. We may specify a frame not only by directly 
specifying vp (the most obvious being vp = 0), selecting some 9-dimensional subspace of A/14, but by specifying any 
variable which can somehow be related to an element of vj? through either Eq. ((144(1 or Eq. 1(142(1 . thus defining that 
element of vf, provided, of course, this results in a non-zero, self consistent solution. It turns out, in fact, that the 
only variables we can't specify in this way are £ a b, Ti. a b, W a b, W a b (the latter being equivalent for these purposes to 
the Zerilli tensor defined later) - these are frame invariant. Below we choose our frame so that 5 elements of v^> are 
zero, thus placing 5 rows of Eq. 1(144(1 into M, and explicitly having our core dynamical system represent the four 
true dynamical degrees of freedom - a four dimensional subspace of A/14. 

To summarise: 

• there are 28 propagation equations for 33 variables, which suggests that there are 33-28=5 'frame' degrees of 
freedom in the choice of the two basis vectors u a and n a ; 

• once we have used the time- independence of the background to harmonically decompose the evolution equations, 
these, combined with the original constraint equations, give rise to a total of 19 linearly independent algebraic 
relations (not all of the 27 evolution and constraint equations - represented by F - are independent); 

• using these algebraic relations to eliminate 19 variables leaves 33 — 19 = 14 variables to be solved for; 

• since the contraint and evolution equations propagate consistently (imposing the constraints 'commutes' with 
the hat derivative) we can be sure that the propagation equations for eliminated variables can be dropped, since 
they will follow from the propagation equations for the remaining variables; 

• we find then that 9 propagation equations (for vd) remain (leaving the 14 — 9 = 5 frame degrees of freedom, v>); 

• finally, since we can choose the 5 frame degrees of freedom more or less arbitrarily we really only have 9 — 5 = 4 
true dynamical propagation equations, as we would expect (to see this, imagine choosing 5 of the 9 elements 
of Vd to be anything at all; then the propagation equations for those variables just fix all five elements of vp, 
so that the only unknowns that remain are the 4 components of V£> for which we have propagation equations. 

A. Determining the full solution: Finding v, M and B 

1. Odd 



a. General Frame If we don't specify a frame choice, and choose our solution vector as, say, V£> = 
(Yij , Qr , Xy , Ay) , then there are two undetermined variables, which we can choose to be vp = (r2 v ,av), as wit- 
nessed in the propagation equation for this solution vector: 

= B g v D + A v F . (145) 



All other variables are linear combinations of elements of v D , except Ev, which depends on Qy, nothing depends on 
av- Thus, Vf represents a frame freedom in the odd-parity variables, and we can specify it at will. 

b. Specific Frame Here, we will choose the frame specifically such that Yy = Ay = 0. This immediately implies 
that £s = = ov = Xy = Zy = fly = 0, which is a rather decent simplification. For our reduced basis vector we 
will choose 



(146) 




i.e., the governing DE will be the Regge- Wheeler equation. The remaining variables are then given by 



V 



Ht 
S t 
Ct 
£v 
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Ay 
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ay 
Xy 
Yy 
Zy 
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—l/iu>4>r 3 















(147) 



where we have used 



C = 3^V - 4 



(148) 



2. Even 



a. General Frame If we don't specify a frame choice, and choose our solution vector as, say, vc E = (Sj, Ct, + 
fly,ay,Xy,Yy), then there are three undetermined variables, which we can choose to be v Fe = (ay , fly, 9s). The 
underlying propagation equation is then 

w D = & g Ey D + A E v F (149) 

All other variables are linear combinations of elements of V£> E , except which depends on Recall that these 
freedoms arise from the nature of the propagation and evolution equations: only the combinations of variables, 
S a — s a b^l b and X — are determined by the propagation and evolution equations, and not the individual variables 
themselves (although they appear in different combinations in other equations); meanwhile, a a is only present on the 
rhs of the propagation equations. Thus, ay, &y and 9$ represent frame freedom in the even-parity variables, and we 
can specify these at will (or other variables which indirectly fix these). 

b. Specific Frame To concur with the odd case above, we will choose here a frame in which Ay = Yy = Zy = 0. 
We choose 




(150) 
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so that the governing DE will be Eq. i|131|) . Thus, the basis vector in V33 can easily be converted to the Regge- 
Wheeler equation, or the Zerilli equation, depending on ones frame of mind [by substituting for Wr or Z from the 
rhs of Eq. (E5) ]. The other variables are linear combinations of the elements of this solution basis vector: 
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(151) 



B. Gravity waves 

The plane wave solutions are given by setting A = 4> = 0, with n a lying in the direction of propagation (although, 
if we keep <fi = 2/r - i.e., set m = - then it represents the same solution, but n a is a radial direction with an 
arbitrary centre). In this case only the four TT tensors are non-zero, and £ ab and Ti. ab represent the curvature of 
plane gravity waves: dynamical tidal forces orthogonal to the direction of propagation. In the BH case the y th erefore 
may be thought of as representing the same part of the GW, although with distortion from the BH itself |24| . This 
is seen in the wave equations which these tensors satisfy which no longer close, and have forcing terms from other 
quantities (in a general frame): 

£\ab } - £{ab] ~ 5 2 £ ab - (0 + 5A) £ {ab} - (0 2 + 8.A 2 - b£) £ ab = 2 (<j> - 2A) 6 {a £ b} + 3£ (0 - 2A) Cab, (152) 

n { ab} - n {ab} - 5 2 n ab - An {ab} - (\<t> 2 - 2A 2 + \£ ) n ab = (y - a) s {a n b} - 3^e c{a s b} c . (153) 

Thus, the principle forcing term for £ ab comes from the shear ofn a , while Ti. ab is principally forced by the shear of 
u a . 

These TT tensors decouple from all the other variables, but not as wave equations, so it is worthwhile giving the 
differential equations which these satisfy. Let W = (Wo, We) where 

*.-(*). and )%-(*). (154) 
Then W satisfies the 4-dimensional first order DE 

/ -v n \ 

(155) 




where 



X = r 2 <P {r 2 (f) 2 [8L - r 2 ((f 2 + 256w 2 )] - 16} , (156) 

-B 2 - 3r 4 cj) 4 (r 2 4> 2 + 4) B x 
-B 1 -ALl{iuy 1 r 2 4> 3 B 2 



T = ( -B 2 - 3r 4 4 (r 2 4> 2 + 4) + 480r 4 </> 2 w 2 -4 (B 1 - 2A0itu 3 r 6 (f> :i ) \ ^ , . 

E ~ v -Si -b 2 J ' 

Si = iwr 2 <?(> [-r 4 (f> 2 (9(t> 2 + 16w 2 ) + 8r 2 2 (2L + 1) - 16] , (159) 
B 2 = -|r 6 6 + 6r 4 2 [<?(> 2 (L - 1) + 64w 2 ] + 8r 2 2 (L - 1) + 32. (160) 

Pretty untidy, but the structure is simple. Note that although decoupling these equations will result in second order 
ODE's for each of the four variables, these are not true wave equations, because of the oj's in T; powers of w up to 3 
are found in Eq. 1155|) . which correspond to third order dot-derivatives. 

In fact, although the Regge- Wheeler tensor, as we have written it, involves £ a fc which is a kinematical term, it may 
be written from purely Weyl contributions, using e.g., Eqs. 1(61(1 and 1)50(1 . or just Eq. 1(51(1 . Therefore, the entire 
solution may be related relatively simply to the Weyl curvature. Indeed, the variables £ a b,'Hab,W a b are all frame 
invariant so the relations between these given for the odd and even cases apply regardless of frame; hence, W a b and 
W a b may be given as linear combinations of E a b and Ti a b, though not in tensorial form. 



The Zerilli tensor 



We have discussed at the end of the section on the Zerilli variable, Sec. lIVBl how it is possible to find an odd parity 
Zerilli variable as a linear combination of Ej, £t and £s, which was seemingly unrelated to the even parity form of the 
Zerilli variable. However, we may utilise our solution for the odd parity variables [Eq. 1)142(1 . although the details of 
this 'general frame' solution are not given here] in order to write Z as a linear combination of Ej, 7iv and 7is, thus 
suggesting a connection with the even parity form, Eq. i(123|) : 

[\L 2 l 2 uj- 2 r^^A- 2 + l]Z = fcg" 1 [3r0E T - 2A~ 1 Hv] - l^^A^Us- (161) 

As yet, the similarity with the even case is not quite clear, although nearly there. If we note that the factor multiplying 
Z on the left hand side is actually a constant (in fact, it is just the special quasi-normal mode factor, discussed earlier), 
then we can re-define our definition of Z, without affecting the differential equation it obeys. Therefore, we define 
the Zerilli tensor Z a b by the differential equation: 

[2 - r 2 (2£ + S 2 )] Z ab = ir0E Qb - §rA- 1 e e{a 8 c H b} - ^S' 2 [2 - r 2 (2£ + 5 2 )] e c{a 5 c S b} U, (162) 

which is related to our odd and even Zerilli variables by 

Z T = [\L 2 l 2 u- 2 r-*r A A- 2 + 1] Z, (163) 
Z T = Z, (164) 

and, of course, both Zy and Zj satisfy the Zerilli equation ((124(1 . 

The complexity of the definition of this Zerilli tensor l(162|) is somewhat startling, especially given the simplicity of 
the analogous Regge- Wheeler tensor. Indeed, the necessity of defining Z ab via a differential equation really does not 
inspire confidence in its possible fundamental nature; were it not for the harmonic splitting we are permitted to do 
due to the nature of the background, we would not be able to actually write down a closed wave equation for Z ab - 
the closed tensorial equation it obeys is 



9 - r 2 <5 2 ) 



Z a b — Z ao — AZ a i 



+ f r - 2 { (4> 2 r 2 - 2 - r 2 8 2 f + 2 (f 2 r 2 - 2 - r 2 ,5 2 ) (|^V - 5 - r 2 S 2 ) 2 } Z ab = 0. (165) 
Hence, only the spherical harmonic amplitudes of Z ao , Zj and Zj, obey a true wave equation. 



D. Discussion 



We have seen in our approach that we may define a TT tensor which satisfies the Regge- Wheeler equation, irre- 
spective of parity, and we have shown that the Zerilli equation may be derived from this tensor in the even case. 
The transformation equations between the two, Eqs. 1(128(1 and 1(130(1 . then allowed us to derive an odd parity Zerilli 
variable, and hence a Zerilli tensor. This transformation between parities is made explicit in our approach due to the 
unification properties of W a b and Z a b- To contrast with the metric approach, the choice of metric functions which 
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makes this similarity between the parities explicit is a physically unmotivated expression, which is a complicated 
linear combination of the metric perturbation functions (see, e.g., Eq. (154) in 0). Of our two fundamental tensors, 
it is fairly clear that the Regge- Wheeler tensor is the most appealing, for two reasons: it is defined in a clear and 
simple way, and it obeys a covariant wave equation. Recall that it's the harmonic amplitudes of Z a b which obey wave 
equations, and not Z a b itself. 

In the frame we have chosen where A a = Y a = Z a = we just have W a b governing both the odd and even 
parity perturbations, which obeys the covariant wave equation H116J) . The Regge- Wheeler equation describes the 
dynamically free gravitational field which propagates at the speed of light. It is given here in its fully covariant, 
gauge- and frame-invariant form. To understand what W a b actually embodies, recall its definition which applies 
in any frame, Eq. (|115Jl . The first term, £ a fc, is just the shear distortion of our sheet (vibrating 2-'spheres') as we 
move radially along n a . The second term is a little more complicated. Recall that in an exact spherically symmetric 
spacetime, £ is the tidal force measured by our static observers u a in the direction n a . There seems to be no reason 
to change this interpretation in the real perturbed spacetime - here it will just undergo fluctuations, but we can still 
call it the radial tidal force. 'Radial' loses its meaning in a perturbed spacetime, but if we chose a frame in which 
n a oc ii a (i.e., A a = 0) an observer could always determine this direction because it would lie precisely in the direction 
of the external force that they must apply. Spatial fluctuations in the radial tidal forces are characterised by (D a £ ) /£ 
which is the comoving fractional gradient of the tidal force; projecting this onto the sheet gives us the gauge invariant 
(first-order) part of this 3- vector: (8 a £) /£ = X a /£, which tells us the fractional gradient of the radial tidal forces over 
a sheet - i.e., how the radial tidal forces change from point to point on our vibrating 2-'spheres'. The distortion of 
this, <5{ a / £\ , appearing in W a b is the shearing distortion of the radial tidal force gradient. The Regge- Wheeler 

tensor is thus of a shearing form; it is this tensor which describes GW around a black hole, through the covariant 
form of the Regge- Wheeler equation (|116|l . 

VI. CONCLUSIONS 

We have presented a new perturbation formalism for dealing with systems with spherical symmetry in the back- 
ground, and we have applied this to the simplest of such systems, the Schwarzschild black hole. Our 1+1+2 splitting 
allowed the Schwarzschild solution to be given in covariant form. We then demonstrated that we can derive the 
usual equations governing perturbations of the spacetime, namely the Regge- Wheeler equation ()117l) . and the Zcrilli 
equation (|124(l . while discussing in detail our new method. We have also shown that there exist Regge- Wheeler and 
Zerilli tensors which unify the odd and even parity perturbations; indeed, the Regge- Wheeler tensor was shown to 
satisfy a closed covariant wave equation which governs the dynamics of the whole problem. This sets the basis for 
future studies of more general astrophysical systems, which only require an appropriate change of the background, 
for which the possible applications are myriad. 

The method has several important aspects. The first is the covariant spacetime splitting itself. In general, the 
two threading vectors u a and n a may be chosen arbitrarily, defining the sheet on which the vectors and tensors 
exist (strictly speaking, in general the sheet is not a true surface, but a collection of tangent planes). This makes the 
approach a halfway house between the 1+3 approach and the orthonormal tetrad approach, and provides a completion 
of the covariant formalism. (Recall that a unique tetrad cannot be defined in isotropic or locally rotationally symmetric 
spacetimes so this is the ideal compromise between the two in such cases.) For systems with spherical symmetry 

in the background, we have seen that when n a is radial in the background, the perturbed spacetime becomes a 
tractable problem, because all vectors and tensors become first order, allowing decomposition with suitable harmonic 
functions - spherical harmonics, in this case. Time harmonics are also introduced to simplify the solution, allowed 
when the background is static, but these are not strictly necessary, as the dot-derivative is a scalar operator, and can 
be dealt with by standard techniques. So far this merely writes the equations in an alternative form. At this point 
finding the solution is relatively trivial: our approach simply requires solving a linear system of algebraic equations, 
and this is all there is to it. The underlying dynamics are then automatically given by the small system of differential 
equations that remain - wave equations, if desired, then may be derived (where they exist) by differentiating this 
equation. 

An important physical aspect of our approach is that it uses a set of covariantly defined quantities with genuine 
physical significance, which makes it clear which objects are crucial for the detection and measurement of gravitational 
waves. Put simply, GW detectors essentially measure gravitational tidal forces; that is, they are sensitive to the 
dynamical Weyl curvature, encoded in the electric Weyl tensor, E a b, and this dynamical Weyl curvature forces the 
working parts of any GW detector through the right hand side of the geodesic deviation equation. We have shown 
that there is a gauge-invariant TT tensor that describes GWs of either parity, and is closely related to the variation 
of the radial tidal force. Thus it is clear that our formalism is dealing with real, physically measurable, objects from 
the start. Indeed, we have discussed how a subset of four of all thirty-three variables is sufficient to determine the full 
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dynamics of the spacetime. 

There are many possible extensions of the work we have presented. The obvious extension is to perturbations of 
static stars, but we envisage that our method will be widely applicable to many other astrophysical situations, such as 
systems with a dynamical background - e.g., collapsing stars, type la supernovae, etc. - where the scalar background 
equations have two (non-tensorial) derivatives in them. In the static background case, where we can introduce time- 
harmonic functions in the perturbed equations, solving the equations is a simple problem of solving a linear system of 
equations, and one is left with a first order non-autonomous system of ordinary differential equations, whose dimension 
is small compared to the original system, plus linear relations among the remaining variables in terms of the basis 
vector of the dynamical system. All the physics of the spacetime is contained in this small dynamical system. With 
a dynamical background, it may not be as simple as this, but we do not envisage it being much more difficult to find 
the full solution; one may have to be careful choosing one's observers (perhaps comoving with the matter in the case 
of a collapsing star). In any event, there is much that can be achieved. In any situation where there is a preferred 
spacelike vector field present, the covariant 1+1+2 sheet formalism should provide new insight. 
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